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ABSTRACT 

We analyze the shapes of the HI velocity profiles of The HI Nearby Galaxy Survey (THINGS) to 
study the phase structure of the neutral interstellar medium (ISM) and its relation to global galaxy 
properties. We use a method analogous to the stacking method sometimes used in high redshift HI 
observations to construct high signal-to- noise (S/N) profiles. We call these high S/N profiles super 
profiles. We analyze and discuss possible systematics that may change the observed shapes of the 
super profiles. After quantifying these effects and selecting a sub-sample of unaffected galaxies, we 
find that the super profiles are best described by a narrow and a broad Gaussian component, which arc 
evidence of the presence of the Cold Neutral Medium (CNM) and the Warm Neutral Medium (WNM). 
The velocity dispersion of the narrow component range from ~3.4 to ~8.6 km s _1 with an average 
of 6.5±1.5 km s _1 , whereas that of the broad component range from ~10.1 to ^24.3 km s _1 with an 
average of I6.8±4.3 km s _1 . Wc find that the super profile parameters correlate with star formation 
indicators such as metallicity, FUV-NUV colors and Ha luminosities. The flux ratio between the 
narrow and broad components tends to be highest for high metallicity, high star formation rate (SFR) 
galaxies. We show that the narrow component identified in the super profiles is associated with the 
presence of star formation, and possibly with molecular hydrogen. 

Subject headings: galaxies: fundamental parameters - galaxies: ISM - galaxies: spiral - galaxies: 
dwarf - ISM: kinematics and dynamics -radio lines: ISM 



1. INTRODUCTION 

Gas in the interstellar medium (ISM) in disk galax- 
ies can exist in a molecular, atomic or ionised phase. 
Most of it will be in the form of neutral atomic hydro- 
gen (HI) whose hyperfine transition can be detec ted at 
the 21 -cm wavelength. Since the early work of |Clark| 
(1965), the neutral atomic ISM has been known to ex- 



alone. 



ist m two phases, known as the Cold Neutral M edium 



(CNM ) and the Warm Neutral Medium (WNM). [Clark 
( |1965[ ) first suggested that HI seen in emission and ab 
sorption had different line widths because it arose from 
two distinct components of the ISM. The narrow absorp- 
tion spectra are caused by a cold component of the ISM 
(the CNM) with a temperature of ~ 100 K whereas the 
broad emission spectra arose from a warm component 
with a temperature as high as 10 4 K (the WNM). The 
existence of these two compo nents was later confirmed by 
Radhakrishnan et al. ( 1972 ) using interferometric obser- 
vations ot HI emission and absorption spectra towards 35 
extragalactic sources situated at intermediate and high 
Galactic latitude. The CNM and the WNM are known 
to coexist in pressure equilibr ium over a narrow range of 
pressure dWolfire et al.||2003[ ). The CNM dominates at 
high densities and pressures whereas the WNM prevails 
at low densities and pressures. These two components 
have been identified in the Milky Way, in nearby dwarfs, 
and spiral galaxies by comparing HI seen in emission and 
absorption as well as from observations of HI emission 
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Young fe Lo| ( |1996 |1997[ ) and |Young et al.| ( |2003[ ) ana- 
lyzed HI emission velocity profiles of seven nearby dwarf 
galaxies. By fitting these profiles with Gaussian compo- 
nents, they found evidence for a broad component with 
a dispersion ranging from about 8 to 13 km s" 1 as well 
as a much narrower component with a dispersion ranging 
from 3 to 5 km s . These velocity dispersions are larger 
than the predicted thermal line widths for the CNM and 
the WNM, showing that other energy sources are also 
playing role in broadening the line profiles. Moreover, 
they found that the narrow component shows a clumpy 
distribution and tends to be located near star forma- 
tion regions, whereas the broad component shows a more 
ubiquitous distribution. Because of the similarities of the 
observed properties of these two components with those 
of the CNM and WNM in M 31 and in our Galaxy, they 
associated the narrow component with the CNM and the 
broad component with the WNM phases of the ISM . A 
similar study was made by |de Blok fc Walter ( 2006 ) for 
the Local Group dwarf galaxy INGG 6822. They also 
found narrow and broad HI components with mean ve- 
locity dispersions of 4 km s _1 and 8 km s" 1 , respectively. 
Here also, the narrow component is usually located near 
star forming regions, whereas the bro ad component tends 
to be found along every line of sight. Braun (19971 ana- 
lyzed HI emission of 11 nearby spiral galaxies and found 
what he called a high brightness filamentary network 
(HBN) of HI emission which he associated with the CNM 
and a diffuse low brightness emission which he identified 
with the WNM. The HBN dominates inside the optical 
radius where it accounts for 60-90 % of the HI emis- 
sion, whereas the diffuse component dominates outside 
this radius. By averaging spectra with peak brightness 
temperature higher than 4er within annular ellipses in his 
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sample, he found that the combined spectra were charac- 
terized by a narrow core superposed on broad Lor entzian 
wings. The width of the narrow line cores allowed Braun 
( 1997|) to giv e an upper limit of 300 K for the HBN. 
Braun (1997) also found that the kinetic t empera t ure of 



the HBJN increases with increasing radius. Braun (19971 
interpreted this as a result of the decrease in the mid- 
plane thermal pressure towards the outer disks. 

Understanding the properties of the gas in the ISM 
is crucial as it is the fuel for star formation. Most of 
this star formation occurs in Giant Molecular Clouds 
(GMCs). These are self gravitati ng clouds wit h masses 
- 10 5 M and sizes ~ 40 pc (e.g. |Tielens|2005[ ). Hydro- 
gen molecules (H 2 ) are the main constituents of GMCs, 
but unfortunately H 2 is not easily observable. Carbon 
monoxide (CO) is the second most abundant molecule in 
GMCs and can be detected at millimeter wavelengths. 
Typically, CO line emission is used to infer the amount 
and distribution of H 2 in galaxies using s ome CO-to-H 2 
conversion factor (e.g., Leroy et al.| |2009). However, the 
strength of the CO emission depends on the meta llicity 




of the ISM flBolatto et "aT]|201l[ |Leroy et aL~|j2011[ ) and 
it is therefore difficult to detect in low metallicity dwarf 
galaxies. An alternative approach is thus necessary to 
investigate the distribution of H 2 there. One example of 
this is to use the dust emission at in frared wavelength s 
as a tracer of the molecular gas (e.g., Leroy et al.||2011 ). 
Another approach, which will be the topic of this paper, 
is to try and identify a precursor of the molecular gas. If 
the H 2 component forms from the cooling of the CNM, 
then identifying the CNM in galaxies can yield indepen- 
dent information on the properties of a component of the 
gas that is closely linked to star formation. 

In this paper, we attempt to detect the CNM and 
WNM components of the ISM by analyzing the shapes 
of stacked HI emission profiles of galaxies observ ed as 
part of The HI Nearby Galaxies Survey (THINGS, [WaT 
ter et al.|[2008l ). THINGS is a large survey of 34 neaTTT 
spiral and dwarf galaxies performed with the NRA 1 
Very Large Array (VLA). The purpose of THINGS was 
to obtain high spectral and spatial resolution maps of 
the neutral hydrogen in these galaxies. The THINGS 
galaxies have distances between 2 and 15 Mpc, star for- 
mation rates of ~ 10 -3 to 10 -6 M Q yr _1 , absolute B- 
luminosities of —11.5 to —21.7 mag, total HI masses of 
(0.01—14) xlO 9 M Q , and metallicities of 7.5 to 9.2 in 
units of [12 + log(0/H)]. Using the VLA B, C and D 
arrays, the THINGS targets were observed at a spatial 
resolution of ~ 6" (robust weighting) or ~ 11" (natural 
weighting), respectively. Two of the THINGS galaxies 
were observed at a spectral resolution of 1.3 km s _1 and 
the rest at a resolution of cither 2.6 or 5.2 km s _1 (see 



Table 2 of |Walter et al. (J2008J)). We use the natural- 
weighted (and non-residual scaled) cubes in this pa- 
per, though our results are independent of the weight- 
ing sc heme used. We refer the reader to |Walter et al.| 
(2008) for an explanation on how the data cubes were 
produced. We use the task PBCORR in the Groningen 
Image Processing System (GIPSY) to correct for the ef- 
fect of the primary beam of the telescope in the image 

4 The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 



cubes. The publicly available THINGS moment maps 
already are corrected for primary beam attenuation. 

This paper is organized as follows. In Section [2 we 
describe the principles of our method for obtaining 
signal-to-noise (S/N) stacked profiles. Section [3 
scribes the fitting and decomposition of the profiles. 
Section |4j we discuss various effects that may influence 
the shapes of the profiles. In Section |5j we describe how 
we define a clean sample based on the studies of these 
effects. In Section|6j we check the reliability of the profile 
shapes of the clean sample. In Section [7J we check pos- 
sible correlations between the shapes ofthe profiles and 
properties of galaxies or morphology. In Section [5J we 
investigate whether the narrow component is associated 
with molecular gas. We summarize our results in Section 

m 

2. GENERATING HIGH SIGNAL-TO-NOISE PROFILES 

Individual velocity profiles in typical HI observations 
generally have a low (< 5) to moderate (< 10) S/N. 
In the THINGS observations, fitting individual, separate 
profiles in a data cube needs to be done carefully to avoid 
skewed results due to the effects of noise. Here, we use 
a method analogous to the stacking meth od sometimes 
used in high-redshift HI observations (e.g., Fabello et al. 
2011 1 to construct high S/N profiles. Note that we stack 
individual spectra within galaxies rather than combining 
global galaxy spectra. Since individual profiles located 
at different positions in a galaxy have different radial 
velocities, the profiles cannot be summed as they are, 
but need to be shifted in velocity to the same reference 
velocity. We use the GIPSY task SHUFFLE to do this. 
This task uses a velocity field to define the amount by 
which a profile in the corresponding data cube needs to 
be shifted. We use Gauss- Hermite polynomial velocity 
fields as input for SHUFFLE. Gauss Hermite polynomial 
velocity fields give robust estimates for the radial velocity 
of the peak of the profile, even in the presenc e of profile 
asymmet ries. For an extensive description see |de Blok et| 
[at| p008l ). 

To minimize the effects of noise, we only include indi- 
vidual profiles whose peak fluxes exceed the 3<r noise level 
in the data cubes. For each galaxy, we sum the shuffled 
profiles to produce what we call super profiles. These 
are high quality profiles, which due to their increased 
S /N allow us to reliably explore the presence of the cold 
and warm neutral components of the ISM. However, as 
information on the individual profiles is lost during the 
stacking, care needs to be taken that systematic effects 
and possible asymmetries in the input profiles are quan- 
tified, to prevent these from mimicking the presence of 
multiple components in the super profiles. This is dis- 
cussed in detail in Section [H 

We show in Figure [T] examples of the super profiles 
of three galaxies (one quiescent galaxy: DD0154, one 
starburst galaxy: NGC 1569, and an interacting galaxy: 
NGC 5194; the rest are shown in Appendix A). These 
super profiles are summed over the entire HI disks of 
the galaxies using the 3cr cut-off mentioned previously. 
Typically, each super profile is the sum of a few hundred 
to a few thousand independent individual profiles (see 
Table [TV. We also plot 3<r uncertainties (as defined in 
SectionpJj), though in many cases these are smaller than 
the symbols used to plot the profiles. The super profiles 
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are characterized by broader wings and narrower peaks 
than purely Gaussian profiles. Most of the super pro- 
files are symmetrical. These properties are all similar to 
what was found for the in dividual profiles of seven ne arby 
dwarf galaxies stud ied by |Young fc Lo (1996 1997]) and 
Young et al. (2003), but observed at a much higher S/N. 

in the following, we will analyze the shapes of the super 
profiles and see if they can be used to infer the presence 
of the CNM and WNM components of the ISM. 

3. PROFILE FITTING AND DECOMPOSITION 



Following de Blok & Walter (2006), we fit the super 
profiles with single Gaussian and double Gaussian com- 
ponents using the GIPSY task XGAUPROF. All param- 
eters (dispersion, central velocity, amplitude, constant 
background term) are left free in the fits. The uncertain- 
ties in the data points of each super profiles are defined 
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where o~ c h.map is the rms noise level in one channel map, 
Nprof is the number of stacked profiles at a certain ve- 
locity V, and N pro f ^ eam is the number of profiles in one 
resolution element (i.e., the number of pixels or profiles 
per beam). The inverse square of these uncertainties is 
used as the weight during the fitting. Note that in Eq. [I] 
we assume that the noise is uniform throughout each data 
cube. Though strictly speaking incorrect (due to the pri- 
mary beam correction), this is a reasonable assumption 
as most of the galaxies only occupy the inner quarter or 
so of the area of the primary beam. The primary beam 
correction factors are therefore on average small. For 4 
galaxies, we correct for the presence of a negative base- 
line level (caused by missing zero-spacings) prior to the 
Gaussian fitting procedure. We fit the baseline with a 
polynomial and substract it from the super profiles. To 
illustrate this method, we show in the Appendix the su- 
per profiles before and after the correction as well as the 
polynomial fit to the baseline. In general the correction 
is of the order of 10% for the velocity dispersions and 
20% for the fitted area of the super profiles. More de- 
tails about this can be found in the Appendix. 

In addition to the one and two component fits, we have 
also tried three- component Gaussian fits to check if this 
further improves the quality of the fits. Comparison of 
the quality of the different fits in Fig. [2] and Table [T] as 
well as the amplitudes of the residuals in Fig. [T] con- 
vincingly shows that the super profiles are non-Gaussian 
and b est desc r ibed by two Gaussian components. Fol- 
lowing Braun (1997), we also fit the super profiles with 
a Lorentzian function. These results are presented in 
the Appendix. Figure [3] shows histograms of the de- 
rived velocity dispersions from both the single and the 
double Gaussian fitting. The distribution of the narrow 
component in Fig. [3] has a mean of 7.6 ± 2.5 km s _1 , 
whereas that of the Droad component has a mean value 
of 22.4 ± 10.0 km s _1 . Also, the broad component has 
a long tail of up to 60 km s -1 . The single Gaussian 
velocity dispersion has a mean of 12.5 ± 3.5 km s _1 . 
If we only consider galaxies with inclinations i < 60°, 
the single Gaussian velocity dispersion has a mean of 
10.9 ± 2.1 km s _1 . This is in agreement wit h the mean 
second-moment value of 11 ±3 km s _1 which |Leroy et al 



( 2008 1 derived for the same galaxies. Table [T] summarizes 
the various fit parameters. 

4. VELOCITY DISPERSION OR BULK MOTIONS 

Even though we have established that the super pro- 
files are best described by two Gaussian components, in- 
terpretation of this is complicated by the possible pres- 
ence of streaming or radial motions and asymmetrical 
input profiles, as well as the limited sensitivity and res- 
olution of the observations. All of this could potentially 
affect the shapes of the super profiles. The brief analysis 
presented in the previous section did not take these ef- 
fects into account. Here we investigate these effects, and 
attempt to identify a sample of galaxies where the su- 
per profiles trace the intrinsic velocity dispersions of the 
ISM. In Appendix B we also present a number of tests 
dealing with the accuracy of the shuffling method itself. 

4.1. Modeling the effects of resolution and inclination 

The gradient in rotation velocity across the finite beam 
of a telescope can broaden velocity profiles. Although the 
THINGS galaxies have been observed at relatively high 
resolution, this could still affect the shapes of the profiles, 
especially for highly inclined galaxies. 

To quantify these effects, we construct model data 
cubes of IC 2574 and NGC 2403 using the GIPSY task 
GALMOD. IC 2574 is a dwarf galaxy with a solid body 
rotation curve, a distance of ^4 Mpc and an inclination 
of 53°. NGC 2403 is a well resolved, rotation dominated 
galaxy typical of spirals found in the THINGS sample. 

For the NGC 2403 models we use the obs erved posi- 
tion angle and inclination (63°) as given in de Blok et 



al 



( |2008 ). We adopt a vertical Gaussian scale height 
of 100 pc, a single component velocity dispersion of 8 
km s , a constant HI surface density and a constant 
rotation curve of 130 km s _1 (the maximum observed 
rotation velocity of NGC 2403). We adopt a channel 
spacing of 5.2 km s _1 as in the THINGS observations. 
We also create an almost identical model, but with an 
inclination of 80° (note that the highest inclination in 
our sample is 76°). We smooth these model data cubes 
to various resolutions and create the corresponding super 
profiles. We fit the artificial super profiles both with a 
single and a double Gaussian to see how the dispersion 
changes with resolution and to check whether resolution 
and inclination can create broad and narrow components 
from a single input component. We have also tested a 
case where the input rotation curve rises linearly up to 
130 km s -1 at 1 kpc and flattens. The results from the 
constant and the linearly rising plus flat rotation curve 
models are identical. 

For the IC 25 74 models, we use t he average position 
angle derived by |de Blok et al. (2008), but assume an in- 
clination of 63° and 80 u to facilitate a direct comparison 
with the NGC 2403 models. Here we also use a single 
component velocity dispersion of 8 km s^ 1 , a constant 
HI surface density but an input rotation curve that rises 
linearly to 80 km s _1 at the edge of the HI disk. We use 
a channel spacing of 2.6 km s _1 as used by THINGS for 
IC 2574. We adopt a vertical Gaussian scale height of 
0.1 kpc. We fit the model super profiles of the IC 2574 
with both single and double Gaussians but the results of 
the double Gaussian decomposition are not meaningful 
as the super profiles are almost perfectly Gaussian. So, 
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Figure 1. Three examples of super profiles for the THINGS galaxies (the rest are shown in Appendix A). For each galaxy, we show the 
results from the single (top panel) and double (bottom panel) Gaussian fitting. The bottom panel of each plot show the residuals from 
the fits. Filled circles indicate the data. The solid lines represent the results from the single and double Gaussian fitting. The dotted and 
the dashed lines represent the narrow and broad components required in the double Gaussian fitting. We plot error bars as 3<r error bars, 
though in most cases they are smaller than the symbols plotted. 
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Figure 2. Solid line histogram: ratio of the reduced x from the 
triple and single Gaussian fits (xIg/Xig)- Dotted line histogram: 
ratio of the reduced x 2 from the double and single Gaussian fits 
(x^q/Xig)' Hatched histogram: ratio of the reduced x 2 from the 
triple and double Gaussian fits (xIg/xIg)' 



Figure 3. Histograms of velocity dispersions derived from both 
one and two-component Gaussian fits. The dotted line histogram 
represents the velocity dispersions derived from the single Gaussian 
fit. The solid line and hatched histograms represent the velocity 
dispersions of the double-Gaussian narrow and broad components, 
respectively. 
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Table 1 

Fitted parameters of the THINGS super profiles 



Galaxy a lg a„ &b deg. A n /A b x\ G lx\c 

asym. 

(km s _1 ) (km s _1 ) (km s _1 ) (km s _1 ) 
1 2 3 4 5 6 7 



DDO 53 


9.9±0.2 


6.2±0.1 


15.0±0.2 


0.2±0.1 


0.46±0 


.02 


0.03 


272 


DDO 154 


9.7±0.1 


6.0±0.1 


13.0±0.2 


0.1±0.1 


0.42±0 


.03 


0.04 


1116 


Ho I 


8.9±0.2 


5.4±0.1 


13.3±0.2 


0.3±0.1 


0.42±0 


.02 


0.04 


380 


Ho II 


9.1±0.2 


5.2±0.1 


12.6±0.2 


0.3±0.1 


0.38±0 


.02 


0.04 


2400 


IC 2574 


10.1±0.2 


5.9±0.2 


14.6±0.3 


0.5±0.2 


0.44±0 


.03 


0.08 


4224 


M81 dwB 


11.3±0.2 


6.6±0.1 


15.2±0.2 


0.2±0.1 


0.28±0 


.01 


0.04 


40 


M81 dwA 


8.2±0.2 


3.4±0.2 


10.4±0.3 


0.1±0.2 


0.14±0. 


.01 


0.12 


88 


NGC 3031 


11.2±0.2 


6.7±0.1 


19.1±0.2 


2.9±0.2 


0.59±0 


.02 


0.05 


8668 


NGC 628 


9.0±0.2 


4.4±0.1 


11.6±0.1 


O.OiO.O 


0.25±0 


.02 


0.02 


7964 


NGC 925 


12.7±0.2 


8.6±0.2 


22.9±0.6 


0.2±0.2 


0.75±0. 


.05 


0.03 


4208 


NGC 1569 


22.5±0.7 


13.0±0. 1 


36.6±1.2 


6.9±0.7 


0.38±0 


.03 


0.09 


440 


NGC 2366 


12.1±0.3 


7.9±0.1 


20.0±0.3 


1.0±0.2 


0.59±0 


.02 


0.03 


1776 


NGC 2403 


11.1±0.3 


6.6±0.1 


18.4±0.2 


0.2±0.1 


0.54±0 


.04 


0.06 


19796 


NGC 2841 


16.3±0.4 


10.4±0.1 


43.2±0.9 


0.8±0.7 


0.72±0 


.03 


0.06 


2600 


NGC 2903 


13.3±0.3 


8.8±0.1 


28.9±0.7 


1.7±0.5 


0.77±0 


.01 


0.06 


2600 


NGC 2976 


11.9±0.3 


8.4±0.1 


20.9±0.4 


1.7±0.3 


0.83±0 


.04 


0.03 


640 


NGC 3077 


12.6±0.3 


5.6±0.2 


17.2±0.2 


0.0±0.2 


0.27±0 


.01 


0.04 


2956 


NGC 3184 


11.3±0.3 


6.1±0.1 


18.8±0.2 


0.6±0.2 


0.43±0 


.01 


0.02 


4860 


NGC 3198 


13.1±0.3 


8.8±0.1 


22.6±0.4 


0.1±0.2 


0.76±0 


.03 


0.03 


2356 


NGC 3351 


10.2±0.3 


7.0±0.1 


22.4±0.6 


0.3±0.4 


0.86±0 


.05 


0.06 


2296 


NGC 3521 


17.4±0.4 


12.4±0.2 


50.6±1.3 


0.4±0.9 


0.89±0 


.01 


0.08 


1448 


NGC 3621 


11.4±0.3 


8.7±0.1 


29.6±1.1 


4.1±0.7 


1.20±0 


.08 


0.08 


4740 


NGC 3627 


20.8±0.5 


14.3±0.2 


43.2±0.8 


0.5±0.5 


0.82±0 


.03 


0.04 


344 


NGC 4214 


8.8±0.1 


4.5±0.0 


11.9±0.1 


0.0±0.1 


0.31±0 


.01 


0.02 


2596 


NGC 4449 


13.8±0.2 


6.9±0.4 


16.8±0.3 


0.4±0.2 


0.21±0. 


.02 


0.04 


2016 


NGC 4736 


11.0±0.3 


7.4±0.1 


22.3±0.6 


0.2±0.4 


0.79±0 


.05 


0.03 


2312 


NGC 4826 


13.1±0.5 


7.9±0.4 


32.0±0.8 


1.6±0.6 


0.61±0. 


.03 


0.04 


264 


NGC 5055 


14.0±0.4 


8.3±0.2 


25.2±0.6 


0.1±0.4 


0.59±0 


.03 


0.06 


6212 


NGC 5194 


17.0±0.4 


9.9±0.1 


29.6±0.4 


1.2±0.2 


0.54±0 


.02 


0.02 


2616 


NGC 5236 


11.2±0.2 


5.6±0.1 


16.8±0.2 


0.7±0.1 


0.38±0 


.03 


0.04 


1912 


NGC 5457 


13.1±0.2 


6.4±0.1 


15.9±0.0 


0.1±0.1 


0.21±0. 


.02 


0.02 


19000 


NGC 6946 


10.4±0.2 


6.1±0.1 


18.3±0.3 


0.4±0.2 


0.55±0 


.03 


0.07 


7252 


NGC 7331 


19.5±0.4 


11.5±0.1 


35.3±0.3 


0.2±0.2 


0.58±0 


.01 


0.01 


2360 


NGC 7793 


10.6±0.2 


6.6±0.1 


17.9±0.2 


0.0±0.1 


0.60±0 


.02 


0.03 


1996 



Note. — Column 1: Name of galaxy; Column 2: Velocity dispersions derived from the one com- 
ponent Gaussian fit; Column 3: Velocity dispersions of the narrow component; Column 4: Velocity 
dispersions of the broad component; Column 5: Degree of asymmetry of the super profiles defined 
as the offset between the peak velocities of the narrow and broad components; Column 6: Ratio of 
the area of the narrow component to the area under the broad component. Column 7: Ratio of the 
X 2 values from the single and double Gaussian fitting; Column 8: Number of independent resolution 
elements used to make super profiles. 



for the IC 2574 models, we only present the single Gaus- 
sian fitting results. We have also tested the case where 
the scale height is 0.5 kpc but the results are similar to 
the 0.1 kpc models. 

Figure [4| and Table [2] summarize the model super pro- 
file parameters as a function of the resolution. It is clear 
that profiles become broader with decreasing resolution, 
and more so for a higher inclination. For example, at 
II" (the average resolution of THINGS), the difference 
between the velocity dispersion measured from the 80° 
and the 63° inclination model data cubes of NGC 2403 
is ~ 9%. At a resolution of 60", the difference is ~ 29%. 
From Figure |4j it is apparent that we do not recover 
the input 8 km s _1 dispersion even at very high resolu- 
tion. This is a result of the finite channel spacing of the 
NGC 2403 data cubes used (5.2 km s _1 ) in combination 
with the high inclination. However, for the IC 2574 mod- 
els, the velocity dispersions approach the input 8 km s _1 
value, and this is due to the smaller channel spacing used 
(2.6 km s" 1 ). 

Figure [5] shows a selection of the modelled super pro- 
files of NGC 2403. It is clear that resolution effects 



can introduce a spurious broad component. We use the 
ratio between the areas of the narrow and broad com- 
ponents A n /Ab (i.e., the flux or mass ratio) to gauge 
whether any of the models resemble the data. Figure [6] 
shows the A n /Ai, ratio as a function of resolution. This 
ratio decreases with resolution, showing the increasing 
importance of the spurious broad component. The ex- 
tremely high ratios at high resolution should not be over- 
interpreted. These are simply due to the fitting routine 
being forced to fit an almost perfectly Gaussian profile 
with two components. At the lowest resolution in the 
models, the broad component created by the low res- 
olution accounts for ~ 13% of the total flux. In the 
THINGS data, the broad component accounts for up to 
~ 70% of the total flux. This discrepancy also indicates 
that resolution effects are not solely responsible for the 
observed broad components. Comparison of the models 
and data shows that at 11" (the typical THINGS res- 
olution), the resulting super profile is well fitted by a 
single Gaussian function. It also shows that the shape of 
the super profile only starts to significantly deviate from 
a single Gaussian at a resolution about six times worse 



G 



14 - 



12 



h 10 




10 20 30 40 50 60 70 
resolution [arcsecl 



13 



12 



w 11 - 



10 




10 20 30 40 50 60 70 
resolution [arcsecl 



26 



24 



" 22 
S 20 
tf 18 




10 20 30 40 50 60 70 
resolution [arcsecl 



Figure 4. Super profile parameters derived from the model data cubes of IC 2574 and NGC 2403 as a function of the resolution. Solid and 
open symbols represent the models of IC 2574 and NGC 2403, respectively. The circle symbols represent the 63° models, whereas the star 
symbols represent the 80° inclination models. We only present the single Gaussian fitting results of the IC 2574 models as these profiles do 
not show a double Gaussian function signature. Left panel: single Gaussian velocity dispersion. Middle panel: Narrow component velocity 
dispersion. Right panel: Broad component velocity dispersion. 

than the THINGS resolution. In summary, models with 
single Gaussian input profiles are well recovered by the 
THINGS resolution. 

The previous test investigated whether double compo- 
nents can be spuriously created using single component 
input models. We now investigate whether input dou- 
ble components can be recovered from limited resolution 
models. We use the 63° inclination model of NGC 2403 
described before. Recall that this model assumes a filled 
HI disk with a Gaussian velocity dispersion of 8 km s _1 . 
We will here refer to this as the broad component model. 
We also construct a second set of cubes, assuming a 
4 km s _1 dispersion and disks that are only partially 
filled (we assume 20%, 40%, 60% and 80% area filling 
factors). We refer to these as the narrow component 
models. We add each of the narrow component models to 
the broad component model, smooth the resulting cubes 
to different resolutions and then derive the corresponding 
super profiles. The narrow component models are scaled 
to give input A n /Ab ratios equal to unity in the regions 
where both narrow and broad components are present. 
Figure [7] shows the measured narrow and the broad com- 
ponents velocity dispersions derived from the smoothed 
models, as a function of resolution. At high resolution, 
the super profile shapes are sensitive to area filling fac- 
tor of the narrow component. As we increase the number 
of narrow profiles in the disk, the super profiles' narrow 
and broad component velocity dispersions approach the 
input dispersions of the narrow and broad component 
disks, respectively. At low resolutions, the super profile 
parameters become insensitive to the narrow component 
filling factor. The models thus show that at the THINGS 
resolution we are able to recover the presence of narrow 
and broad components from the models. 



4.2. Effect of high intensity profiles 

Another possible explanation for the non-Gaussianity 
could be that the shapes of the super profiles are dom- 
inated by a few lines of sight with high intensities and 
narrow dispersions. This implies some correlation be- 
tween peak brightness of the input profiles and velocity 
dispersion. To test this, we create super profiles using 
the 50% brightest and 50% faintest (in terms of peak 



brightness) profiles in each galaxy. 

We compare the parameters of the faint and the bright 
super profiles of each galaxy in Figure [8j Most veloc- 
ity dispersion points scatter around the fine of equality. 
However, four high dispersion galaxies systematically de- 
viate. These are NGC 2841, NGC 3521, NGC 3627 and 
NGC 925 (in order from highest to lowest difference). 
This may indicate some kind of bulk motion of gas caused 
by, e.g., interaction or tidal effects in these galaxies. In 
terms of area, the faint super profiles show A n j ratios 
that tend to be smaller than those of the bright super 
profiles. The mean difference between the bright and 
faint A n /Ab ratios is ~36%. As most of the faint profiles 
are found in the outer disks of the galaxies, this seems 
to indicate that the A n /Ai, ratio in the outer HI disks is 
smaller than in the inner disks. This result will be dis- 
cussed further in a future paper. On the whole though, 
it is clear that the dispersions of the super profiles of the 
bright and faint parts are not significantly different. 

As an additional test, we also create super profiles by 
normalizing individual profiles by their peak intensity be- 
fore summing them. This way, any possible effects of a 
small number of bright input profiles are removed. Fig- 
ure [8] shows that, in terms of velocity dispersion, the 
normalized super profiles are very similar to the original 
ones. The correlation coefficient is R 2 ~ 0.99 for both 
the narrow and the broad components. In terms of area, 
the A n /Ai, ratios derived from the normalized super pro- 
files tend to be smaller than those from the original super 
profiles. The difference between the mean values is, how- 
ever, only ~1.3%. In summary, the shapes of the super 
profiles are not dominated by high intensity, narrow dis- 
persion input profiles. 

4.3. Effects of radial motions 

Having ruled out resolution effects and high-intensity 
input profiles as the cause for the non-Gaussianity of the 
observed super profiles, we now investigate the effects of 
large-scale radial motions. A net radial motion in one 
direction will create an asymmetric tail in the velocity 
profiles. More specifically, for a given radial motion, the 
location of the tail with respect to the profile depends 
on the spatial position of the profile: a radial motion 
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Table 2 

Effect of the resolution on the shapes of the super profiles 



NGC 2403 models IC 2574 models 





Resolution 




Velocity 




Flux ratio 


Resolution 


Velocity 








dispersion 




A n /A b 




dispersion 
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ClG 


(") 


(pc) 


(km s^ 1 ) 


(km s _1 ) 


(km s ) 




(PC) 


(km s -1 ) 


Inclination: 63° 


4 


62 


8.6 


8.5 


15.5 


83.3 


78 


8.1 


8 


124 


8.6 


8.6 


15.4 


60.6 


155 


8.1 


11 


171 


8.6 


8.6 


15.4 


56.6 


213 


8.1 


20 


310 


8.9 


8.8 


15.6 


23.4 


388 


8.2 


30 


465 


9.3 


9.0 


16.1 


12.4 


582 


8.3 


40 


621 


9.8 


9.3 


16.8 


7.9 


776 


8.1 


50 


776 


10.1 


9.6 


17.5 


5.6 
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8.6 


60 


931 


10.6 


9.8 


18.1 


1.2 


1164 
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Inclination: 80° 
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16.6 


6.6 
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9.6 


17.7 


4.0 
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8.2 


30 


465 


11.4 


10.0 


19.0 


2.6 


582 


8.3 


40 


621 


12.5 


10.6 


20.5 


1.8 


776 


8.5 


50 


776 


13.6 


11.2 


21.9 


1.4 


970 


8.7 


60 


931 


14.9 


11.8 


23.3 


1.2 


1164 


9.0 



Note. — &ig'- Velocity dispersion derived from one-component Gaussian fit. a n : Velocity disper- 
sion of the narrow component, erf,: Velocity dispersion of the broad component. A n /Ab\ Ratio of the 
fluxes or masses associated with the narrow and bro ad c omponents. Here, only the single Gaussian 
results of the IC 2574 models arc shown (see Section I47TI. 
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Figure 5. Observed and simulated super profiles of NGC 2403 derived after smoothing the model data cubes to the resolution indicated 
on the plots. The left and right panels represent the results from the model data cubes of NGC 2403 with input inclination of 63° and 80°, 
respectively. The dotted and the dashed lines represent the narrow and broad components required in the double Gaussian fitting. The 
solid lines represent the results of the double Gaussian fitting. 
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Figure 6. Flux ratio of the narrow and broad components, 
A n /Ai,, derived from the model data cubes of NGC 2403 as a func- 
tion of the resolution. The open circle and star symbols represent 
the 63° and the 80° inclination models, respectively. 
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Figure 7. Velocity dispersions d erive d from the modelled narrow 
and broad disk models (see Sect. |4.1| l as a function of the resolu- 
tion. Open circles represent the narrow component. Filled circles 
represent the broad component. The percentage indicates the nar- 
row profile disk filling factor. 



causing a tail on the high- velocity side of a profile that is 
located on one side of the major axis of the galaxy, will 
create a tail on the low-velocity side of the profile when 
it is located on the opp osite side of the major axis. See 
Warner et ai] ( |1973[ ) and |van der Kruit fc Allen] ( | 19 78 1 ) for 
a detailed discussion of these kind of patterns typically 
observed in velocity fields and profiles. 

Figure [9] shows examples of super profiles extracted 
from the opposite sides of the major axes of each galaxy 
overplotted on top of the overall super profiles of these 
galaxies. We use a consistent definition of sides corre- 
sponding to their average position angle (PA). One side 



of the galaxy is chosen as the one with an average posi- 
tion angle of 90° with respect to the receding semi-major 
axis (denoted as the PA+90° side); the other side has an 
average PA of 270° with respect to the receding semi- 
major axis (the PA+270 side). We show in Figure 10 
a histogram of the difference between the velocity cfTsT 
persion derived from the opposite sides with respect to 
the major axes of each galaxy. For about 90% of the 
studied galaxies, the difference is less than 1.5 km s _1 . 
There are four galaxies with higher velocity dispersion 
differences. These are NGC 3031, NGC 3077 (both inter- 
acting), NGC 4826 (counter-rotating HI disk) and NGC 
7331. The latter has the highest inclination (76°) in the 
sample and the difference in velocity dispersion for this 
galaxy may be a projection effect. The velocity disper- 
sion comparison presented above only works if the dis- 
tortion due to radial motions is asymmetric. Therefore, 
we also measure the skewness of each of the super pro- 
files on the opposite sides of the major axis. We fit the 
super profiles with a Gauss-Hcrmitc h% function. The 
skewness coefficient is defined as 4 v / 3/i3 and the results 
are shown in Fig. |11| The apparent trend visible in the 
figure is mostly caused by highly inclined, interacting or 
disturbed, as well as starburst galaxies. It is possible 
that these galaxies suffer from radial motions, and thus 
they will not be part of our subsequent analysis, as de- 
scribed in detail in Sect.[5j The remaining outliers have 
very symmetrical profiles with only some faint extended 
wings showing up at levels well below 20% of the peak 
intensity. While these affect the skewness, they contain 
negligible flux, and have no impact on the dispersions or 
flux ratios. Most likely these apparent wings are caused 
by residual baseline effects. 

In the following, we develop simple models to under- 
stand how radial motions affect the shapes of the super 
profiles. We do this by creating models which have dif- 
ferent ratios of rotation and radial motions. We again 
use t he 6 3° inclination NGC 2403 model discussed in 
Sect. 4.1 We create four new model data cubes by 
adding outward radial motions of 5, 10, 15, and 20 
km s _1 respectively to the original model (thus shifting 
the input profiles in velocity) . We again assume a single- 
component Gaussian velocity dispersion of 8 km s , a 
constant rotation curve and a constant surface density 
HI disk for all models. 

For each of the 5, 10, 15 and 20 km s _1 model data 
cubes, we randomly choose 20%, 40%, 60%, 80% of the 
profiles and add these to the original NGC 2403 model 
(which contains no radial motions). We thus create 16 
data cubes with different fractions of the disk area con- 
taining asymmetric profiles with different amounts of ra- 
dial motion. We smooth these model data cubes to the 
THINGS 11" resolution. The velocity fields derived from 
the combined, smoothed cubes are then used to shuf- 
fle profiles and derive the corresponding super profiles, 
which we th en fit with single Gaussian functions. 

Figure 12 shows the resulting velocity dispersions (as 
measured with a single Gaussian) as a function of the 
input expansion velocities. As expected, the velocity 
dispersions increase with increasing expansion velocities. 
The increase is steeper as we increase the fraction of gas 
with non-zero expansion velocities. We see that a ra- 
dial motion of 5 km s -1 causes only minimal broaden- 
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Table 3 

Effect of the radial motions on the 
shapes of the super profiles 



Vrad 






Mb 




(km s~ 1 ) 


20% 


40% 


60% 


80% 


5 
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99.3 
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12.9 
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Figure 8. Top panel: comparison of super profile parameters derived from the faint and bright parts of galaxies. Bottom panel: comparison 
of super profile parameters from the normalized and the original super profiles. The triangle shaped symbols represent interacting or 
disturbed galaxies (NGC 3031, NGC 4449, NGC 5194, NGC 5457, NGC 3077). The open circle symbols represent galaxies with star bursts 
or mild AGN activity (NGC 1569, NGC 3521, NGC 3627). The star symbols indicate non-interacting galaxies that have an anomalously 
high velocity dispersion (NGC 7331, NGC 2841). The rest of the galaxies are represented by the filled circle symbols. 

ing, independent of the area coverage. The most ex- 
treme model with 20 km s _1 radial motions and 80% 
area coverage has a superprofit that is broadened by 
~ 23%. Large-scale non-circular motions of this mag- 
nitude are rare in disk galaxies, and mostly only found 
in the central parts of barred galaxies. We can com pare 
this with the work by Trachternach et al. (2008) who 
used harmonic decompositions of the velocity fields of 
the THINGS galaxies to quantify the non-circular mo- 
tions there. They found a median amplitude of all non- 
circular motions of A r = 6.7 ± 5.9 km s _1 . Of their 18 
analyzed galaxies, 17 have A r < 10 km s _1 and only one 
has A r > 20 km s _1 . 

Finally, we repeat the procedure described above, but 
rather than assuming the same constant HI surface den- 
sity disk for all models, we assume that the profiles 
with non-circular motion have an amplitude that is 30% 
of that of the profiles without non-circular motions. 
This choice will allow the non-circular motion profiles 
to mimic low-level broad wings in the super profiles. In 
this case we fit all resulting s upe r profiles with double 
Gaussians. We show in Figure [12] and Table [3] the fitted 
An/A}, ratios as a function of input expansion velocities. 
The A n /Ab ratio decreases with increasing expansion ve- 
locities. If we assume that the non-circular motio ns in 
the THINGS galaxies are of order 10 km s _1 or less (Tra- 



Note. — V ra d- Input radial velocity used 
to construct the model data cubes. Percent- 
age indicates the area coverage of gas with 
non-circular motions. A n /Ab'. Flux ratio of 
the narrow and broad components. 



chternach et al. 2008), then these motions are only able 
to introduce spurious broad wings containing at most 7% 
of the total flux, i.e., much smaller than observed. We 
therefore conclude that radial motions cannot be a ma- 
jor cause of the observed non-Gaussianity of the observed 
super profiles. 

4.4. Thick disks and lagging halos 

Lastly, we investigate the possibility that our sample 
galaxies have a thin disk embedded in a "lagging" thick 
HI disk with a lower rotation velocity, as was found for 



10 



Typical 



3 



-a 
B 



0.14 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 

j 

I 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 

NGC 2976 J 




i i 1 i 
NGC 7793 J 


i i i i 1 i i i i 1 


i i i i 1 i i i i 
NGC 2366 J 

- 

\ 


III II II 1 


III II I I 1 






1 1 1 1 1 1 1 1 


Mill 






-1— 1 — |— 1 — 1— 1 — h 


-| — 1— 1 — I I 






" " 




" 


" " 



-150 -100 -50 50 100 150 -100 100 -100 -50 50 100 



Offset velocity [km s '] 



Extreme 



X3 
CD 
N 



O 



0.5 




0.14 


-0.14 



1111,1 


111,1111,1111 




, , | , i , , , , , 


111,1111, 






NGC 4826 J 




NGC 3077 J 




NGC 3031 J 








1 V 




f \ 












III 


1 1 1 1 1 1 1 1 1 1 1 


-1— 1 — ^-r— 


+— 1 — 1 ,.„l — I I 

, , \"~ 


1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


, , , 1 '",'< , , 1 , , , , " 




, , ■ , 1 , , , , 1 , , , , 1 , , ■ 



-200 -100 



100 200 -100 100 -150 -100 -50 50 100 150 

Offset velocity [km s _1 ] 



Figure 9. Super profiles extracted from the opposite sides of the major axes of each galaxy. The super profiles from the two halves of 
each galaxy are represented as dashed and dotted lines. The overall super profiles are represented by the solid lines. At the bottom of each 
panel is the difference between the super profiles from the two sides and the total super profiles. Note that the super profiles are normalized 
to their peak values. 

have differences less than 1 km s ; for the remaining 
10 galaxies, the difference is between 1 and 4.5 km s _1 . 
Of those 10 galaxies, 4 are interacting (NGC 3031, NGC 
4449, NGC 5194, NGC 5457), 2 are kinematically dis- 
turbed (NGC 1569 and NGC 3521), one is known to be 
a galaxy with a counter-rotating inner disk (NGC 4826). 
To complement our analysis, we also measure the coeffi- 
cient of skewness of the super profiles in the approaching 
and receding sides. The result is p resented in Figure 
[l5[ Most of the points in Fig. [15] are offsets towards 
the lower right quadrant. This means that super pro- 
files in the receding sides tend to be negatively skewed, 
whereas those in the approaching sides tend to be posi- 
tively skewed. These are consistent with but do not by 
themselves prove the existence of lagging thick HI disks. 
Most of the outliers in Fig. [15] will not be considered for 
further analysis as describedm Sect. [5] 



NGC 891 (Sw atgrsetalJlJ^llOosterloo et al.|2007[ ) and 
NGC 2403 flFraternali et al.||2002[) . The super profiles 
then represent a "cross-cut" through the global dynamics 
of the system, rather than a reflection of the state of the 
ISM. If this were the case, the lagging disk would broaden 
the profiles asymmetrically depending on whether the 
approaching or the receding side of the galaxy is studied. 
Adding these two halves would give the impression of the 
super pr ofile having broad wings. 

Figure [13] shows examples of super profiles from the 
approaching and receding sides of a few typical sample 
galaxies plotted on top of the super profiles derived from 
the entire disk. Note that the super profiles plotted there 
are normalized to their peak values. For reference we also 
show the most asymmetric approaching a nd r eceding su- 
per profiles in the bottom panel of Fig. [13] We com- 
pare the velocity dispersions derived from the approach- 
ing and receding sides in Figure [14] For each galaxy 
the largest of the velocity dispersions derived from the 
two halves of each galaxy is shown as ai arge , whereas the 
smallest one is given by a sma ii ■ We also show a histogram 
of the difference between the vel ocity dispersions derived 
from the two sides in Fig. [14] Of the 34 galaxies, 24 



5. DEFINING A CLEAN SAMPLE 

In this section we define a subsample of galaxies that do 
not suffer from the effects of interaction or major bulk 
motions. Figure 16 shows the measured super profile 
velocity dispersions plotted against the inclinations of 
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Figure 10. Histogram of the difference between the velocity dis- 
persion derived from the opposite sides of the major axis, where 
PA+90° side refers to the position corresponding to a position an- 
gle (PA) of 90° with respect to the PA of the receding major axis 
and PA+270 side is the position corresponding to a PA offset of 
270° from the receding major axis. 
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Figure 11. Skewness (defined as 4\/3/i3, where /13 is the third 
coefficient of the Gauss-Hermite function.) of the super profiles 
derived from the opposite sides of the major axis. PA +90° side: 
position corresponding to a position angle (PA) of 90° with respect 
to the PA of the receding major axis. PA+270 side: position cor- 
responding to a PA offset of 270° from the receding major axis. The 
three outlying galaxies are, from left to right, NGC 1569, DDO 53 
and NGC 3031. The plus and minus signs represent the signs of 
the skewness parameters. Symbols follow Fig. [8] 

the THINGS galaxies. We see a narrow relation with 
galaxies with higher inclinations h aving broader pro files 
(as also shown for the se galaxies in Leroy et al.|2008| and 
Tamburro et al.|2009 1 . Many of the galaxies that do not 
follow the relation (NGC 1569, NGC 5194, NGC 2841, 
NGC 3627, NGC 3521, and NGC 7331) are known to be 
affected by strong sources of turbulence or bulk motion 
(such as interaction with other galaxies or starbursts). 



These galaxies are marked with different symbols in Fig. 
16 One thing to note is the remarkable constancy of the 
Cn/ffb ratio. The weighted average value of this ratio is 
0.37 ±0.04. 

Based on inspection of Figures [10} [TT| fl4| |l~5| and [To] we 

exclude from the rest of our analysis galaxies, that ao not 
follow the observed t rend between inclination and veloc- 
ity dispersion in Fig. [l6j We also exclude galaxies where 
the difference in velocity dispersion measured for the ap- 
proaching and receding sides, or PA+90° and PA+270 
sides is more than 1.5 km s . Interacting galaxies will 
also not be considered. Based on these criteria, we are 
left with 22 galaxies which we identify as our clean sam- 
ple. Their super profiles are expected to be less domi- 
nated by the various effects presented earlier. The clean 
sample is listed in Table |4j 

6. EFFECTS OF ASYMMETRIC PROFILES 

Having identified a sample likely not to be suffering 
from large-scale disturbances, we now take a detailed 
look at the shapes of the observed input profiles in this 
sample. 

If the individual profiles used in the stacking are them- 
selves asymmetric, this can also lead to the observed su- 
per profiles, as the sum of a negatively skewed and a 
positively skewed profile will be one with broader wings 
and a narrower peak than a Gaussian profile. 

Here we check whether the broad components in the 
super profiles are intrinsic or artifacts due to asymmetric 
input profiles. We do this by selecting only symmetrical 
profiles, creating super profiles from those symmetrical 
input profiles and checking whether broad and narrow 
components are still present. 

To select symmetrical profiles, we create masks based 
on the difference between the third order Hermite 
(HER3) velocity field values and intensity weighted mean 
first moment values (IWM). The former trace the veloc- 
ity of the peak of the profile, the latter the intensity 
weighted mean value. For asymmetric profiles there will 
therefore be a difference between the HER3 and the IWM 
values. For symmetrical profiles the difference between 
HER3 and IWM is small, and this can be used to effi- 
ciently select symmetrical input profiles. If super profiles 
created from only the symmetric profiles still show broad 
components, then these are likely to be intrinsic. We thus 
define three masks based on the difference in velocity be- 
tween HER3 and IWM, A H -i, as follows: 

— Symmetrical profiles (SP): positions where |A#_/| < 
5 km s . 

— Left-handed asymmetrical profiles (LHAP): positions 
where A H _j < —5 km s -1 . 

— Right-handed asymmetrical profiles (RHAP): posi- 
tions where Ajj-i > 5 km s . 

SP therefore selects symmetrical profiles; the other two 
masks LHAP and RHAP select only asymmetrical pro- 
files. We show for all galaxies the location of the sym- 
metrical and asymmetrical profiles in Fig. [T7J We create 
super profiles using the three masks and fit them with 
both single and double Gaussian components. Examples 
of the super profiles of NGC 3521 created using the three 
kinds of masks, overplotted on top of its total super pro- 
files are shown in Fig. [TSl It is clear that the masks 
produce the expected results. 

We compare the dispersions from the original super 
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Figure 12. Left panel: single Gaussian velocity dispersions derived fro m th e first set of model data cubes of NGC 2403 as a function 
of the input expansion velocities used to construct the models (see Sect. [473^ . Right panel: A n /A^ ratio derived from the second set of 
model data cubes of NGC 2403 as a function of the expansion velocities (see sect, |4.3[ l. Percentage indicates the fraction of profiles having 
non-zero expansion velocities in the model data cubes. 

profiles (i.e., those derived from using the entire HI disk; 
see Sect. [2| with those derived using the SP mask in 
Fig. 19 (top panel). The dispersions of the narrow com- 
ponent from the original super profiles are similar to 
those from the SP super profiles. However, the disper- 



sions of the broad component from the original super 
profiles tend to be larger than those from the SP super 
profiles. This is what we would expect if the original pro- 
files were broadened by asymmetry in the input profiles. 
We can remove this effect by only selecting symmetric 
profiles. 

To test whether the |A#_r| < 5 km s" 1 criterion pro- 
duces symmetrical profiles, we create new masks with 
|A#_j| < 2km s _1 . This selects only very symmetrical 
profiles but at the cost of a reduced number of input pro- 
files. Using the more strict condition results in a loss of 
about 60% of the profiles. Figure 19 (bottom panel) com- 
pares the |A#_/| < 5km s~ x SP mask velocity disper- 
sions with the |A_ff_/| < 2km s _1 equivalent. It is clear 
that the dispersions of the broad components of the lat- 
ter tend to be smaller by ~ 1 km s _1 . Despite the smaller 
num ber of input profiles in the |Ajj-j| < 2 km s _1 case, 
Fig. [19] shows that the uncertainties have not increased 
dramatically. In the following, we will therefore use the 
2 km s _1 condition unless stated otherwise. 

Results from the single and double Gaussian fitting 
are summarized in Table [4l This convincingly shows 
that both narrow and broad components are needed to 
fit the profiles. We also fit the wingless sides of LHAP 
and RHAP super profiles (i.e., the "Gaussian" half of 
the profile with single and double Gaussian components. 
Here we also find that the two-components fit perform 
better. 

The broad component velocity dispersions derived 
from the wingless side of the LHAP and RHAP super 
profiles tend to be smaller than those derived from the 
SP super profiles. However, there is no major difference 
between the narrow component velocity dispersions de- 
rived from the wingless side of LHAP and RHAP super 



profiles and those from SP super profiles (see Fig. 20). 
The measurements of the narrow component velocity dis- 
persion are thus more robust than those of the broad 
component. 

7. SUPER PROFILES AND GLOBAL PROPERTIES OF 
GALAXIES 



Figure 21 shows the distribution of the narrow and 
broad component velocity dispersions derived using the 
SP profiles in the clean sample galaxies. The narrow 
component has a much tighter distribution than the 
broad component, with mean values of 6.5 ± 1.5 km s _1 
and 16.8 ± 4.3 km s _1 , respectively. We adopt the defi - 
nition of spirals and dwarfs given in |Leroy et al. ( 2008 ) 



They define dwarfs as galaxies with rotation velocities 
Vrot < 125 km s _1 , stellar masses M* < 10 10 M , and 
absolute B magnitude Mb > —20 mag. Galaxies more 
massive than this are defined to be spirals. Dividing 
our clean sample into dwarfs and spirals, we find their 
velocity dispersions are identical within the uncertain- 
ties. For the dwarfs, we find a mean narrow compo- 
nent velocity dispersion of 6.5 ±1.5 km s _1 and a mean 
broad component velocity dispersion of 15.1±3.6 km s . 
For the spirals, we find values of 6.7 ± 1.4 km s" 1 and 
18.7 ± 4.3 km s _1 , respectively. 

The A n /Af, ratios differ though. For the spirals we find 
a mean of 0.63±0.24; for the dwarfs a value of 0.49±0.20, 
indicating a larger importance of the narrow component 
in spirals. We summarize the fitted parameters of the 
clean sample super profiles in Table [4] 

In the following, we explore possible correlation be- 
tween global galaxy properties and super profile shapes. 
If the narrow and broad components identified here are 
related to the CNM and the WNM phases of the ISM, 
we expect some correlation between these two compo- 
nents and tracers of star formation. For example, |Schaye| 
( 2004 1 suggests that the transition from the warm to the 
cold phase of the ISM triggers disk instabilities which 
eventually leads to star formation. We use FUV-NUV 
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Figure 13. Super profiles extracted from the approaching and receding sides of the galaxies. The super profiles from the two halves of 
each galaxy are represented as dashed and dotted lines. The overall super profiles are represented by the solid lines. The bottom panels 
of each plots represent the difference between the super profiles from the two sides and the total super profiles. The super profiles are 
normalized to their peak values. 

colors, Ho luminosities and metallicities as star forma- 
tion indicators. The FUV-NUV colors are drawn from 



Lee et al.|(2011[ ) or else |Muaoz-Mateos et al.| ( |2009p . The 
(FUV-INUV) colors have been corrected lor galactic ex- 
tinction. Emission in the UV arises from O and later- 
type B stars and therefore traces star formation over a 
timescale o f < 100-200 Myr s . The Ha luminosities are 
taken from Kennicut et al. (2008). Emission in Ha is 



mainly caused by hot stars such as O and early-type B 
stars which are massive and short lived. Thus, Ha emis- 
sion is a tracer of recent star formation over a timecale 
of < 10 Myrs. This timescale is roughly c omparable to 
the ti me it takes for the WNM to cool (e.g. I Wolfire et al.| 
2003|) . The metal abundance is taken from Moustakasl 
et al. l pOTol) or alternatively from Table 1 of Walter et 
P005f~[Moustakas et at] (|2010[) used both empirical 



is likely to be underestimated by up to ~0.2-0.3 dex in 
the empirical calibrations and overestimated by a simi- 
lar amount in the theoretical calibrations. Thus, taking 
the average values of the two calibrations is expected to 
yield an abundance estimate that is close to the "true" 
value. Met allicity plays an important role in cold phase 
formatio n ( |Schaye||2004| |Krumholz et al.||2009| |Walch et 
al. 2011) as gas can be cooled, tor example, via collisional 



ilyugin & Thuan 2005 ) and theoretical ( Kobulnicky & 



|Kewley||2004p calibrations to compute gas-phase a bun 
dances. Here we use the average values of the charac- 
teristic (globally averaged) m etallicities of the two cal - 
ibrations listed in Table 9 of M oustak as et al. (2010) 



As discussed in Moustakas et aT. (2010), the abundance 



excitation of heavy elements such as oxygen. 

The relationship between the super profile parameters 
and these different star formati on indicators are pre- 
sented below. We plot in Figure [22] the velocity disper- 
sion ratio (cr„/<7;,) and the A n /A}, ratio of the narrow and 
broad components as a function of metallicities, (FUV- 
NUV) colors and Ha luminosities. 

We find that the velocity dispersion ratio decreases 
with increasing metallicity with a correlation coefficient 
R ~ —0.56. In principle this could reflect more efficient 
cooling, but as turbulence must also play a large role 
in determining the value of the velocity dispersion, it is 
more likely that this trend is a reflection of a metallicity- 
star formation rate relation. High metallicity galaxies 
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Figure 14. Top left panel: Comparison of the widths of the super profiles from the approaching and receding sides of the galaxies. The 
largest of each of these velocity dispersions is represented by uiarge an d the smallest one by 0" S mall- The dashed line represents the line of 
equality. Top right panel : Histogram of the difference between the velocity dispersion derived from the approaching and receding halves of 
th e galaxies. Bottom left pan el: Narrow component velocity dispersion as a function of the amplitudes of non-circular motions, A, derived 
by Trachternac h et al. 1 2008 . Bottom righ t panel: Broad component velocity dispersion as a function of A. Note that only the 18 galaxies 



that are part ot | Irachternaeh et al.| | |2008[ l sample are shown here. Symbols are as in Fig. [8] 

tend to have higher star formation rates, and for these 
high star formation rate (high metallicity) galaxies, the 
broad component velocity dispersion is likely to be larger 
which will give a lower a n / crj . 

The a n ja h ratio is higher for bluer (FUV-NUV) col- 
ors; the correlation coefficient is i? ~ — 0.61. As low 
metallicity galaxies are on average also bluer in (FUV- 
NUV), this relation probably reveals similar information 
as the metallicity vs <7 n j <j\> relation. There is only a weak 
correlation between dispersion ratio and Ha luminosities 
(B, ~ —0.20). This weak correlation may be due to the 
fact that the super profile parameters and the Ha lu- 
minosities have not been derived from the same regions 
within the galaxies. This will be an interesting topic for 
follow-up studies. 

In terms of the A n /Ab ratio, which should measure the 



amount of cold gas (CNM) relative to warm gas (WNM), 
we find that the fraction of the narrow components rel- 
ative to the broad components tends to be high for high 
metallicity galaxies (R ~ 0.53). We find similar relations 
in terms of (FUV-NUV) colors and Ha luminosities with 
correlation coefficients R ~ 0.53 and R ~ 0.47, respec- 
tively. These correlations are all preliminary indications 
that the presence of cold gas is associated with star for- 
mation. 



S. 



DOES THE NARROW COMPONENT TRACE 
MOLECULAR GAS? 



Star formation and molecular hydrogen are intimately 
connected. Unfortunately, molecular hydrogen is not di- 
rectly observable and CO emission is usually used as a 
tracer. Usually a constant CO-to-H2 conversion factor is 
assumed to derive the amount of H2 present in a galaxy. 
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Table 4 

Fitted parameters of the clean sample 



Galaxy 



1 



Tig 

(kms- 1 ) 
2 



(kms 1 ) 
3 



(kms -1 ) 
4 



An/Ab xl G /xlc 



0.41±0.03 


< 


0.01 


0.50±0.04 


< 


0.01 


0.39±0.02 




0.02 


0.33±0.02 


< 


0.01 


0.42±0.03 




0.10 


0.11±0.01 




0.01 


0.27±0.02 




0.03 


0.74±0.06 




0.02 


0.58±0.03 


< 


0.01 


0.60±0.06 




0.02 


0.78±0.04 


< 


0.01 


0.86±0.05 


< 


0.01 


0.44±0.01 


< 


0.01 


0.64±0.03 




0.07 


0.95±0.06 




0.02 


0.71±0.06 




0.16 


0.29±0.01 




0.25 


0.99±0.07 




0.03 


0.55±0.05 




0.17 


0.36±0.03 




0.03 


0.65±0.09 




0.02 


0.63±0.02 




0.21 



DDO 53 
DDO 154 
Ho I 
Ho II 
IC 2574 
M81 dwA 
NGC 628 
NGC 925 
NGC 2366 
NGC 2403 
NGC 2903 
NGC 2976 
NGC 3184 
NGC 3198 
NGC 3351 
NGC 3621 
NGC 4214 
NGC 4736 
NGC 5055 
NGC 5236 
NGC 6946 
NGC 7793 



9.7±0.1 
9.6±0.1 
8.9±0.1 
8.9±0.1 
9.8±0.2 
8.2±0.1 
8.8±0.2 
12.6±0.2 
12.2±0.2 
10.7±0.2 
12.6±0.3 
11.9±0.2 
10.4±0.3 
13.1±0.2 
9.9±0.2 
11.3±0.2 
8.3±0.1 
10.0±0.2 
13.3±0.3 
10.6±0.2 
10.1±0.2 
10.4±0.2 



6.0±0.1 
6.3±0.1 
5.4±0.1 
5.0±0.1 
5.7±0.2 
3.4±0.1 
4.5±0.1 
8.6±0.2 
8.2±0.1 
6.6±0.1 
8.2±0.2 
8.5±0.1 
5.9±0.1 
8.5±0.1 
7.0±0.1 
7.5±0.2 
4.2±0.0 
7.1±0.1 
7.9±0.3 
5.5±0.1 
6.4±0.3 
6.7±0.1 



13.3±0.2 
12.9±0.2 
12.2±0.2 
11.7±0.2 
13.7±0.3 
10.1±0.2 
11.6±0.1 
20.4±0.8 
17.7±0.3 
16.8±0.2 
24.3±0.7 
19.9±0.4 
17.9±0.2 
20.0±0.3 
21.5±0.7 
18.7±0.7 
10.9±0.1 
20.3±0.7 
20.8±0.7 
15.1±0.1 
17.5±0.1 
17.0±0.2 



Note. — Column 1: Name of galaxy; Column 2: Velocity dispersions 
derived from the single Gaussian fit; Column 3: Velocity dispersions of the 
narrow component; Column 4: Velocity dispersions of the broad component; 
Column 5: A n /At ratio. Column 6: Ratio of the x 2 values from the single 
and double Gaussian fitting. Spiral galaxies (adopting the definition of |Leroy| 
|et al.|[2008| l are marked in bold face; the rest are dwarf galaxies (using the 
same definition). 
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Figure 15. Skewness (defined as 4\/3/i3, where /13 is the third 
Gauss-Hermite polynomial coemcient) of the super profiles derived 
from the receding and approaching side of each galaxy. The two 
outlying galaxies in the lower left quadrant are DDO 53 (top) and 
NGC 1569 (bottom). The plus and minus signs represent the signs 
of the skewness parameters. Symbols are as in Fig. [8] 

However, the strength of the CO emission line depends 
on the metal abundance of galaxies and CO is difficult 
to detect in low metallicity environments. CO may also 
not trace all the H2 content of galaxies because unlike 



dances where both CO and H2 are present, UV photons 
from young stars will easily dissociate CO while leaving 
H 2 intact. 

As H2 forms as the cold phase of the neutral medium 
cools, we test here whether the narrow component identi- 
fied in the super profiles can be used to infer the amount 
of H2 in galax ies. This is motivated by ear ly observa- 
tionsmadebv |Young fe~Lo| ( |1996[ |1997| and fde Blok~fc 
Walter (2006). They found that the narrow component 
tends to be located near regions of star formation. The 
narrow component could thus be associated with the star 
fomation in a manner similar to that of H2. Also, from 
a theoretical point of view, HI must pass through a cold 
phase before turning molecular. 

We consider the 11 galaxies in the clean sample that 
have also been detected in CO in t he HERA C O-Line 
Extragalactic Survey (HERACLES, |Leroy et al 



2009) 



H 2 , CO cannot self-shield (e.g., Krumholz et al. 20111 



and therefore, in regions with low metal and dust a bun 



HERACLES provides maps of CO J = 2 — > 1 emission in 
18 THINGS galaxies (though the full HERACLES sam- 
ple is larger) at 13" angular resolution and 2.6 km s _1 
velocity resolution using observations with the IRAM 30 
m telescope. The CO-detected galaxies have metallicities 
ranging from 8.25 to 9.12 in units of [12 + log(0/H)]. 

If the narrow component is indeed associated with 
molecular gas then we expect correlations between the 
A n /Ab ratio and the H 2 to HI mass ratio. When masses 
are measured in identical regions, this translates in cor- 
relations between the respective surface densities as well. 
We also test for correlations with the narrow and broad 
components velocity dispersion ratios. 

We derive super profiles and average HI and H 2 mass 
surface densities in the overlap regions in the sample 
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Figure 16. Relation between inclination and measured velocity dispersion. All error bars are 3 a error bars. The top left panel is for 
the narrow component. The top right panel is for the broad component and the bottom left panel is for the single Gaussian component. 
The dashed lines are linear least squares fit with slopes of 0.045±0.014 for the single Gaussian component, 0.057±0.010 for the narrow 
component and 0.18±0.05 for the broad component. Note that only the filled data points were used in these fits. The bottom right panel 
is a plot of the ratio of the narrow and broad component velocity dispersions against inclination. The dashed line represents the weighted 
average value of a n /a b ~ 0.37±0.04. Symbols follow Fig. [8] 



galaxies where both CO and HI are detected. Restrict- 
ing ourselves to the overlap regions means we are testing 
only the regime where both gas phases are present. We 
do not investigate the regime where all gas has turned 
molecular and no HI is present anymore (such as the cen- 
tral parts of some galaxies where central HI deficiencies 
are filled with molecular gas). 
The values for Em and Eh 2 are calculated using 



correct for the presence of helium. In the above Equa- 
tion, i?2i is a CO J = 2— s-ltoJ=l— s-0 line ratio 



Em 
Sh 2 



4.4 



X, 



(2) 

cos i / C o(2^1) (3) 

where Em and Ejj 2 are in units of Mq pc~ 2 , and /hi and 
Ico are both in K km s _1 . The mass surface densities 
are corrected for inclination and include a factor 1.36 to 



0.02 cosi /hi 
X 1 

co R21 



which we assume to be equal to 0.7 ( Schruba et al.|2012 ), 
X is a conversion factor normalized to the Galactic CO 
( J = 1 — > 0)-to-H 2 conversion factor X co = 2.0 x 10 20 
cm~ 2 (K km s -1 ) -1 . Here we adopt X = A"co- 

In Figure [23] we compare \og(A n /A b ) and log(cr n /ab) 
with the values for log(E H2 /E H i) as derived for the same 
areas in each galaxy. The observed upper limit on a n /ab 
and A n /Ab seems to vary systematically with Eh 2 /Ehi- 
The respective correlation coefficients are R — 0.47, 
and R = —0.70. A possible interpretation of the ob- 
served trend between \og(a n /a b ) and log(E H2 /E H i) is 
that galaxies having larger fractions of molecular gas are 
expected to have more abundant star formation, which 
drives turbulence to the ISM and results in higher u b val- 
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Figure 17. Location of symmetric and asymmetric profiles for the THINGS galaxies. Blue pixels represent symmetric profiles (SP). Green 
and red pixels represent left-handed asymmetric profiles (LHAP) and right-handed asymmetric profiles (RHAP), respectively. Ellipses 
represent the optical radius r25. The small grey circles indicate the approaching sides. 




Figure |17| (continued). 




1 7 (continued) . 
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Figure 18. Examples of super profiles of NGC 3521 using LHAP (left panel), SP (middle panel) and RHAP (right panel) represented by 
the dashed lines, overplotted on top of the total super profile (shown as solid lines). 
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Figure 19. Top panels: Comparison of the velocity dispersions of the original super profiles to those derived from SP profiles using 
|A^_/| < 5 km s" 1 . Bottom panel: Comparison of the velocity dispersions derived from SP profiles using the |Aj^_/| < 5 km s _1 with 
those derived from SP profiles using |A^_/| < 2 km s _1 . Dashed lines are lines of equality. Symbols are as in Fig. [5] 
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Figure 20. Comparison of the velocity dispersions derived from the SP super profiles to those derived from the wingless sides of the 
LHAP (toppanel) and RHAP (bottom panel) super profiles using |A£f-i| < 2 km s — The dashed lines are lines of equality. Symbols are 
as in Fig. [8| 



ues and thus lowers the \og(a n /o~b) ratio. The trend in 
A n /A\, ratios might reflect the fact that atomic gas passes 
though the CNM phase before turning into molecular. 
Thus, galaxies having larger A n /Ab (more cold atomic 
gas) ratio are expected to have larger molecular gas frac- 
tions. This result serves as early indication that the nar- 
row component may be associated with molecular gas. 
It will be interesting to study the spatial distribution of 
these two gas components in more detail and investigate 
whether the locations of the narrow components corre- 
late with those of the molecular gas. 

9. SUMMARY 

We have conducted in-depth analyses of the shapes of 
the HI velocity profiles of galaxies from The HI Nearby 
Galaxy Survey (THINGS) sample. To minimize the ef- 
fect of noise on the velocity profiles, we have constructed 
high signal-to-noise (S/N) profiles by aligning individual 
profiles in velocity and stacking them. We call these su- 



per profiles. We have quantified the relevant systematic 
effects that may change the intrinsic shapes of the super 
profiles and defined a clean sample where the observed 
shapes of the super profiles are less affected by these ef- 
fects. 

We have fitted the super profiles with single, double 
and triple Gaussian models, as well as Lorentzian mod- 
els. Based on a \ 2 analysis, we found that the shapes of 
the super profiles are optimally described by the sum of 
a narrow and a broad Gaussian component. The narrow 
component velocity dispersions of the clean sample range 
from 3.4 to 8.6 km s _1 with a mean of 6.5 ± 1.5 km s _1 . 
The broad component velocity dispersions range from 
10.1 to 24.3 km s" 1 with a mean of 16.8 ± 4.3 km s _1 . 
Note that due to the limitation of the data, the derived 
velocity dispersions are overestimated by at most 20%. 
The combined effects of the finite channel spacing, the 
limited spatial resolution (see Figure [4]) and the incli- 
nation (see Figure 16) contribute about 15% while the 
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Figure 21. Histograms of the velocity dispersion of the clean 
sample. The solid line histogram represents the narrow component. 
The dotted line histogram represents the broad component. 



effects of radial motions (whose magnit ude is expected 
to be less than ~10 km s _1 , see Section 4.3) contribute 
about 5%. 



We have investigated possible correlations between the 
shapes of the super profiles and global properties of 
galaxies. We found that the flux ratio of the narrow and 
broad components tends to be high for high metallicity, 
high star formation rate galaxies. The flux ratio also 
increases with increasing bluer FUV-NUV colors. We 
interprete these correlations as evidence that cold HI gas 
is associated with star formation. In addition, the veloc- 
ity dispersion ratio of the narrow and broad components 
decreases with increasing metallicity, FUV-NUV colors 
and Ha luminosities. 

We present tentative evidence that the narrow com- 
ponent is associated with molecular gas. We find that 
upper limits on a n ja\, and A n /Ab change with Eik/Shi, 
suggesting that the cold neutral phase and the molec- 
ular phase are related. The observed trend between 
the A n /Ab ratio and the Eh^/Ehi reflects the fact that 
atomic gas passes through the CNM phase before turn- 
ing into molecular. Based on this preliminary analysis, 
it is expected that the location of the narrow compo- 
nent (CNM) correlates with the location of molecular 
gas. Higher resolution pixel-by-pixel studies of the HI 
line profiles should be able to confirm this. 
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Figure 22. Top panel: velocity dispersion ratio as a function of metallicity, (FUV-NUV) colors and Ho luminosities. Bottom panel: Flux 
ratio of the narrow and broad components as a function of metallicity, (FUV-NUV) colors and Ha luminosities. The triangle symbols 
represent interacting galaxies (those that are tidally interacting with a nearby companion or are being affected by their environment; these 
are NGC 3031, NGC 4449, NGC 5194, NGC 5457, NGC 3077). The open circle symbols represent galaxies that are disturbed kinematically 
due to the effect of e.g., star formation (NGC 1569, NGC 3521, NGC 3627). The star symbols indicate non-interacting galaxies that have an 
anomalously high velocity dispersion (NGC 7331, NGC 2841). The diamond shaped symbols represent non-interacting and non-disturbed 
galaxies where the dispersions measured on the opposite sides of the galaxies differ by more than 1.5 kms -1 . (NGC 4826 and M81 DwB; 
see Section [5] for the choice of this value). Filled circle symbols represent our clean sample galaxies. 
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Figure 23. Velocity dispersion (left panel) and A n /Af, ratio (right panel) as a function of H2 and HI mass surface density ratio. 
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APPENDIX 

A. SUPER PROFILES: GAUSSIAN AND LORENTZIAN 

In this appendix, we present the super profiles of the entire THINGS sample (see Section [2] for an explanation on 
how these profiles were generated). Note that we correct for the presence of a negative bowim the super profiles of 
NGC 628, NGC 2403, NGC 5236 and NGC 5457 by fitting a polynomial to the baseline and substract the polynomial 
from the super profiles. Also note that NGC 628, NGC 2403, and NGC 5236 are among the clean sample galaxies (see 
Section [5]). In general, the negative bowl correction is small for these galaxies (of the order of 10% for the velocity 
dispersions and 20% for the fitted area) and their broad component parameters do not appear as outliers when plotted 
against other quantities such as inclination (see Figure 24 for illustration). We summarise in Tableland Table [| the 



super profil es' parameters before and after the negative bowl correction 



Following Braun (1997), we have also fitted the super profiles with a Lorentzian function. We compare the reduced 
X 2 values from the double and Lorentzian fitting in Table[7j Of the 34 fitted super profiles, 4 can be well fitted by either 
a double Gaussian or a Lorentzian function (i.e., their reduced x 2 values agree within 10%). For 25 super profiles, 
a double Gaussian function is clearly preferred. Finally, for 5 super profiles, a Lorentzian function best describe the 
shapes of the profiles. Overall, a double Gaussian function seems to be an optimal description of the profiles. We do 
not attempt to make any physical interpretation of the Lorentzian fitting results, however we summarize the fitted 
parameters in Table [7 

The super profiles of the THINGS galaxies are shown in Fig. |25j For each galaxy, the left panel represents a plot of 
the super profile fitted with a single Gaussian component whereas the middle and right panels represent super profiles 
fitted with a double Gaussian and a Lorentzian functions, respectively. The filled circles indicate the data points. 
The solid black lines represent the results from fitting with a single and double Gaussian, as well as Lorentzians. The 
dashed and the dotted lines represent the narrow and broad components required in the double Gaussian fitting. We 
show the residual from the three kinds of fits at the bottom panel of each figure. For the four galaxies with negative 
bowls mentioned earlier, the open and the filled circle symbols represent the data points before and after the correction, 
respectively. The dashed lines represent the polynomial fits to the negative bowls. 
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Table 5 

Velocity dispersion of super profiles before and after the 
negative bowl correction 



GALAXY 


°1G °Vi 






(kms -1 ) (kms -1 ) 


(kms -1 ) 




Before negative bowl correction 




NGC 628 


8.4 4.0 


10.0 


NGC 2403 


10.4 6.0 


15.0 


NGC 5236 


10.7 5.1 


15.0 


After negative bowl correction 


NGC 628 


9.0 4.4 


11.6 


NGC 2403 


11.1 6.6 


18.4 


NGC 5236 


11.2 5.6 


16.8 



Note. — <?ig'- Velocity dispersion derived from the single 
component Gaussian fit. a n : Velocity dispersion of the nar- 
row component. <Jb'- Velocity dispersion of the broad compo- 
nent. 



Table 6 

Comparison between the fitted area of super profiles before and 
after the negative bowl correction 



GALAXY 


A\G,b.c./ A lG,a.c. 


-^n,b.c. 1 ' -^-n.a.c. 


Ab : b.c. /^-b,a.c. 


NGC 628 


0.89 


0.67 


0.93 


NGC 2403 


0.91 


0.76 


0.94 


NGC 5236 


0.95 


0.81 


0.97 



Note. — b. c. — Before correction (i.e. before negative bowl cor- 
rection), a. c — After correction (i.e. after correcting for the negative 
bowl). 
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Figure 24. Broad component velocity dispersion as a function of inclination. The solid symbols represent our clean sample. The grey 
symbols represent the velocity dispersions of (from left to right) NGC 628, NGC 5236 and NGC 2403 before and after the negative bowl 
correction. 
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Table 7 

Parameters of Lorentzian fit to the THINGS super profiles 



Galaxy FWHM Area mean xIg/Xl. 

(km s _1 ) (Jy km s _1 ) (km s _1 ) 
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Note. — Column 1: Name of galaxy; Column 2: Full Width at 
half maximum ; Column 3: Integrated area; Column 4: Fitted central 
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Gaussian and Lorentzian fitting. 
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Figure 25. Super profiles of the THINGS galaxies. For each galaxy, we show the results from the single (left panel) and double (middle 
panel) Gaussian fitting. The right panels show a fit with a Lorentzian function. The bottom panels of each set of plots show the residuals 
from the fits. The filled circles indicate the data. The solid black lines represent the results from the single, the sum of the double Gaussian, 
and the Lorentzian fitting. The dotted and the dashed lines in the middle panel represent the narrow and broad components required in 
the double Gaussian fitting. We plot error bars as 3a error bars, though in most cases they are smaller than the symbols plotted. 
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Figure |25| (continued). 
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Figure |25| (continued). 
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Figure |25| (continued). 
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Figure |25| (continued). 
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Figure |25| (continued). 
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Figure [25] (continued). Here the open and the filled circle symbols represent the data points before and after the baseline 
correction, respectively (see Section |3j. The dashed lines represent the polynomial fits to the negative bowls. 
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Figure |25| (continued). 
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B. SYSTEMATIC SHUFFLE UNCERTAINTIES 

Choosing the wrong offset velocity in the profile shuffling procedure could artificially broaden super profiles. As the 
SHUFFLE method as implemented in GIPSY uses velocity field values as input, the type and accuracy of the velocity 
field used is important as well. To get an idea of the uncertainties in velocity field values in real data, we show in Figure 
[26| the uncertainties in the centr al velocity value of t he fitted third-order Hermite /13 profiles of NGC 2403 plotted 



against the fitted amplitudes (see de Blok et al.| [2008). As expected, low amplitude profiles have larger uncertainties. 




log(Amplitude) [mjy beanT 



Figure 26. The uncertainties in the central velocity value as derived by fitting a Hermite (/13) function to the profiles in NGC 2403 
plotted against the fitted amplitude of these profiles. 



These uncertainties could in principle create an "artificial" broad component in the super profiles. We therefore carry 
out an experiment to quantify the effect of these uncertainties on the super profiles. We test a case where all profiles 
have random offsets with respect to their proper "shuffle" values as well as a case where the offsets depend on the 
amplitude of the spectra. To do this, we create artificial data cubes. The data cubes contain purely Gaussian profiles, 
which all have 6 km s" 1 velocity dispersion but with different amplitudes. We use the observed peak flux distribution 
of NGC 2403 as input probability distribution for our models. We then randomly pick one thousand amplitudes using 
this distribution and generate Gaussians using these amplitudes and the 6 km s _1 constant velocity dispersion. The 
position along the velocity axis of the model profiles is taken from the observed velocities of the profiles of NGC 2403. 
We create a data cube using the input spectra described above and derive a velocity field from it, which we refer to 
as the true velocity field. We then add uniform random offsets to the true velocity field and use the modified velocity 
field to shuffle profiles. 

For the first test, we give all the input spectra uniform random offsets (between -5 to 5 km s _1 , -10 to 10 km s _1 , 
etc.). For the second model test, only profiles with amplitudes less than 25% of the overall peak amplitude are given 
uniform random offsets (for the case of NGC 2403, where the maximum profile peak value is 9.2 mjy, this value 
correspo nds to 2.3 mjy). Larger amplitude profiles are assumed to have zero offset. 

Figure 27 (left panel) shows examples of super profiles derived from the first set of model cubes. The super profiles 
get broader as we increase the offsets but the two components retrieved by the double Gaussian fitting are all similar in 
amplitude and dispersion. In other words, there are no broad and narrow components. The results of a single Gaussian 
fit to the simulated super profiles show that at [-5, 5] km s _1 offsets, the resulting broadening due to incorrect shuffling 
is about 1 km s _1 . At [-10, 10] km s _1 , this increase to ~ 3 km s _1 . At [-15, 15] km s _1 offsets, the width of the 
resulting super profiles is twice as large as the width of the input profiles. At [-20, 20] km s _1 offsets, the super profiles 
start to have double peaked features. 

Figure 27 (right panel) shows examples of super profiles derived from the second model cube. We start to clearly 



see a broad component resulting from incorrect shuffling at [-15, 15] km s _1 offsets. The wings of the simulated super 
profiles get more and more pronounced with increasing offsets. 
We summarize in Table [8] the velocity dispersion of the super profiles derived from the two tests. 
The above results show that we need an offset of at least 15 km s _1 to explain the non-Gaussianity of super profiles 



by incorrect shuffling. As Figure 26 shows, the maximum uncertainties in the fitted amplitudes are around 5 km s 



This is three times smaller than the offset required to create a broad component. At [-5, 5] km s" 1 offsets, i.e., equal 
to the maximum uncertainty in the real data, the broadest of the two components required in the double Gaussian 
fitting has a velocity dispersion of ~ 6.7 km s _1 at these offsets. Thus, we introduce a broadening of only ~ 10% to 
the broad component by going from to [-5, 5] km s -1 offsets. Obviously, as this is the maximum uncertainty in the 
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Table 8 

Velocity dispersions derived from simulated super profiles using different offsets 



Offsets uniform amplitude dependent 



0~1G a„ tJb °~1G °"n a b 

(km s _1 ) (km s _1 ) (km s _1 ) (km s _1 ) (km s _1 ) (km s _1 ) (km s -1 ) 



-5, 5 7.1 6.5 6.7 6.6 6.1 6.7 

-10, 10 9.2 7.1 7.2 6.9 6.2 7.6 

-15, 15 11.9 8.0 7.8 7.3 6.1 9.4 

-20, 20 15.0 8.8 9.3 7.3 6.1 11.4 

-25,25 ••• ■•• ••• 7.3 6.2 14.6 

-30,30 ••• ••• ••■ 7.2 6.2 18.5 



Note. — (Column 1): Range of offsets. Column 2-4'- Derived dispersions using the 
offset range given in Col. 1 as applied to all profiles. Column 5-7: Derived dispersions 
using the offset range given in Col. 1 as applied only to profiles with an amplitude less 
than 2.3 mjy. <Tig denotes the velocity dispersion derived from a single Gaussian fitting; 
a n the velocity dispersion of the narrow component and (75 the velocity dispersion of the 
broad component. 

data, the real broadening is likely to be much smaller. 

Uncertainties in the SHUFFLE procedure are not sufficient to create the broad component that we see in our super 
profiles. The offsets required to create a broad component are much larger than inferred for real data. 




Figure 27. Examples of simulated super profiles derived by giving all input spectra (left panel) and low amplitude spectra (right panel) 
uniform random offsets. The dotted and the dashed lines represent the narrow and broad components required in the double Gaussian 
fitting. The solid black lines represent the results from the double Gaussian fitting. 



